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EXECUTIVE  SUMMARY 


THE  MAGNETIC  INDUCTION  OF  THE  SYSTEM  CONSISTING 
OF  A COIL  AND  A FERROMAGNETIC  SPHERICAL  BODY 


OBJECTIVE 

The  objective  of  this  work  was  to  derive  solutions  to  static  ferro- 
magnetic problems  that  included  both  current-carrying  coils  and  linear 
ferromagnetic  bodies.  The  solutions  are  intended  for  comparison  with 
solutions  to  ferromagnetic  problems  obtained  by  various  numerical 
techniques  such  as  the  finite  difference  method,  the  finite  element  method, 
and  the  integral  equation  iterative  solution  method. 

APPROACH 

After  deriving  the  governing  differential  equation  from  Maxwell's 
equations  for  classical  magnetostatic  field  theory,  the  method  of  separation 
of  variables  was  employed  to  obtain  the  problem  solution. 

RESULTS 

The  magnetic  induction  was  calculated  for  two  geometries  (configura- 
tions) of  a ferromagnetic  spherical  body  surrounded  by  a current-carrying 
conductor.  The  first  case  was  for  an  infinitesimally  thin  current  band 
carrying  a stationary  current  and  surrounding  a spherical  shell.  The 
second  case  was  for  a current  band  of  finite  width  carrying  a stationary 
current  and  surrounding  a solid  sphere.  The  ferromagnetic  bodies  were 
assumed  to  be  linear  and  homogeneous.  The  reduction  of  the  solution  of  the 
dipole  term  to  that  of  a filamentary,  circular,  current-carrying  conductor 
is  shown  for  the  first  case  when  the  permeability  of  the  ferromagnetic 
spherical  shell  is  allowed  to  approach  that  of  free  space. 


RECOMMENDATIONS 

It  is  recommended  that  the  derived  solutions  be  programmed  on  a 
digital  computer  for  direct  comparison  of  these  results  to  those  obtained 
by  various  numerical  methods.  There  are  plans  to  implement  these  recom- 
mendations during  the  fiscal  years  1979  and  1980. 
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The  magnetic  induction  is  calculated  for  two 
configurations  of  a ferromagnetic  spherical  body 
surrounded  by  a current-carrying  conductor.  The 
first  case  is  for  an  infinitesimally  thin  current 
band  carrying  a stationary  current  and  surround- 
ing a spherical  shell.  The  second  case  is  for  a 
current  band  of  finite  width  carrying  a stationary 
current  and  surrounding  a solid  sphere.  The  ferro- 
magnetic bodies  are  assumed  to  be  linear  and  homo- 
geneous. The  reduction  of  the  solution  for  the 
dipole  term  to  that  of  a filamentary  circular 
current-carrying  conductor  is  shown  in  the  first 
case  when  the  permeability  of  the  ferromagnetic 
spherical  shell  is  allowed  to  approach  that  of 
free  space. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  Program  Element  1121N,  Project  B0005, 

Task  Area  B0005-SL-001,  Work  Unit  2704-110.  The  project  director  is 
Mr.  J.  L.  Corder,  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center. 
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INTRODUCTION 

In  the  past,  exact  analytical  solutions  of  Maxwell's  equations  using 
classical  formulations  have  been  limited  to  body  shapes  and  inhomogeneities 
that  conform  to  a few  separable  coordinate  systems.  With  the  advent  of 
modern  digital  computers  with  large  computational  and  storage  capabilities, 
many  electromagnetic  field  problems  of  importance  in  engineering  practice 
can  and  have  been  solved  by  using  a numerical  solution  to  the  governing 
differential  or  integral  equations  under  a suitable  choice  of  boundary 
conditions.  Such  numerical  solutions  of  Maxwell's  equations,  when  used 
with  a complete  description  of  the  electric  and  magnetic  sources  and  the 
constitutive  laws  of  the  media,  can  be  used  to  describe  completely  the 
electric  and  magnetic  fields  produced  by  the  source,  including  nonsymmetric 
geometries,  nonsymmetric  source  distributions,  and  spatially  varying 
media  parameters. 
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The  literature  contains  numerous  examples  that  demonstrate  the 
power  of  using  numerical  techniques  such  as  the  finite  difference  method, 
the  finite  element  method,  and  the  integral  equation  iterative  solution 
method.  Both  two-  and  three-dimensional  electrical  engineering  problems 
have  been  solved.  Authors  have  tried  to  validate  the  numerical  methods 
and  calculations  by  comparing  numerical  results  to: 

1.  Simple  problems  using  engineering  approximations. 

2.  Other  numerical  method  solutions  of  the  same  problem. 

3.  Laboratory  experimental  results. 

4.  Problems  that  have  analytical  solutions. 

Another  technique  that  could  be  used  but  which  the  authors  of  this  report 
have  not  seen  very  widely  used  in  the  literature  is  a comparison  of  two 
theoretical  formulations  for  a given  problem  solved  by  the  same  numerical 
method.  An  example  might  be  a vector  potential  formulation  compared  to 
a scalar  potential  formulation. 

The  motivation  for  this  work  arose  out  of  the  need  for  solutions  to 
static  ferromagnetic  problems  that  could  be  used  for  comparison  with 
numerical  methods.  The  capability  exists  to  calculate  the  magnetic  field 
due  to  a ferromagnetic  body  immersed  in  a constant  inducing  field  for  the 
following  bodies:  solid  sphere,  spherical  shell,  solid  prolate  spheroid, 
prolate  spheroidal  shell,  solid  general  ellipsoid,  and  general  ellipsoidal 
shell. ^ Also,  the  capability  exists  to  calculate  the  magnetic  field  caused 

g 

by  a single  or  combination  of  current-carrying  coils.  However,  the 
capability  did  not  exist  to  calculate  the  magnetic  field  caused  by  a ferro- 
magnetic body  in  the  presence  of  the  field  created  by  a current-carrying  coil. 


BASIC  EQUATIONS 

We  can  start  with  Maxwell's  equations  for  classical  magnetostatic 
field  problems 


1 


A complete  listing  of  references  appears  on  page  63. 
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T 


Vx  H = J (la)* 

7 ' t = 0 (lb) 

where  H magnetic  field  intensity  (A/m)** 

B magnetic  flux  density  (T  or  Wb/m  ) 

2 

J electric  current  density  (A/m  ) 

In  the  general  case  for  ferromagnetic  materials  B is  a nonlinear  function 
of  it 


B = f (H) 


(2) 


where  as  shown  in  Figure  la  B is  not  a single  valued  function  of  H.  The 
function  f(H)  depends  on  the  magnetic  history  of  the  material,  that  is 
how  the  metal  attained  its  magnetization.  This  is  referred  to  as  hysteresis. 
It  is  also  noted  that  any  property  of  a ferromagnetic  material  has  meaning 
only  if  it  is  considered  together  with  its  complete  magnetic  history. 

In  certain  practical  engineering  problems,  the  variation  in  the 
magnetic  intensity  is  small,  and  the  functional  relationship  between  B and 
H is  approximately  linear  (see  Figure  lb) . For  the  linear  case  where  the 
material  is  isotropic,  the  magnetic  induction  B is  related  to  the  field 
intensity  H by  the  relationship 


3 = UQ  (xj-l)  H = MoMr  H = yH 

where  Xm  magnetic  susceptibility  (dimensionless) 

y magnetic  permeability  (henry/meter) 
(X^+l)  = y^  relative  permeability  (dimensionless) 


yQ  free  space  permeability  (4ttx10 


.-7 


(3) 


henry/meter) 


*The  del  operator  V is  defined  with  respect  to  the  rectangular  coordinate 
system  and  is  strictly  valid  in  a rectangular  coordinate  system  only. 

Very  often y*  and  are  used  as  equivalent  symbols  for  curl  and  divergence 
generally.  TTiis  use  is  followed  in  this  report. 

**Definitions  of  abbreviations  are  on  page  v. 
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(b)  Curve  for  a Ferromagnetic  Material 
at  Low  Inducing  Fields 
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FLUX  OCNUTY I.  (mim1) 


This  report  assumes  that  the  ferromagnetic  body  has  isotropic  and 
linear  material  properties.  The  divergenceless  nature  of  the  magnetic  flux 
density  in  conjunction  with  the  fact  that  the  divergence  of  the  curl  of 
any  vector  function  is  zero  allows  the  introduction  of  the  magnetic 
vector  potential  field  (A) 


^ 

B = V X A (4) 

where  A is  the  magnetostatic  vector  potential  function  in  weber /meter. 

The  substitution  of  equation  (4)  into  equation  (la)  gives  the  fundamental 
equation  of  the  vector  potential  of  the  magnetostatic  field. 

- V x (VxA)  - (VxA)  xyJ=J  (5) 


For  homogeneous  materials  as  assumed  in  this  report  the  magnetic  permeability 
is  spatially  invariant.  Hence 

V J * 0 (6) 

and  equation  (5)  reduces  to 

V x V x A = pJ  (7) 


Using  the  vector  identity 

yxyxt=  V (V-T)  - (8) 

equation  (7)  becomes 

V (V-A)  - $A  = n J (9) 

The  magnetostatic  vector . potential  is  characterized  by  the  important 
property  that  its  divergence  can  be  conveniently  chosen  to  be  zero. 

V.T  = 0 (10) 

Equation  (9)  reduces  to  the  vector  Poisson's  differential  equation. 

£ A = - pj"  x (11) 


This  is  the  governing  equation  for  our  calculations. 


The  general  boundary  conditions  to  be  satisfield  at  the  interfaces 
of  stationary  dissimilar  media  may  be  derived  from  the  limiting  integral 
forms  of  Maxwell's  equations  and  are  given  by 

• (^-B^)  * 0 or  Bnl  = Bn2  (12a) 


x (H2-H1) 


Js  or  Ht2"Htl 


J 

s 


(12b) 


where  the  subscripts  1 and  2 indicate  the  media  under  consideration,  and 
"n^2  denotes  the  unit  normal  vector  to  the  interface  and  is  directed  from 
medium  1 into  medium  2.  In  the  case  where  the  materials  are  linear  and 
isotropic  equations  (12a)  and  (12b)  become 


(12c) 

(12d) 


Jg  is  a true  surface  current  density  that  may  exist  at  the  interface.  At 
an  interface  where  Jg  is  0,  equations  (12b)  and  (12d)  need  to  be  modified 
accordingly. 


THIN  COIL  SURROUNDING  A FERROMAGNETIC  SPHERICAL  SHELL 


GENERAL  SOLUTION 

We  now  proceed  to  solve  the  boundary  value  problem  of  a ferromagnetic 
spherical  shell  of  outer  radius  Rj,  Inner  radius  R^,  and  a homogeneous 
permeability  y_,  surrounded  by  an  infinitesimally  thin  current  band  of 
radius  R^  having  a current  density  J.  A constant  current  density  is 
assumed.  Figure  2 identifies  the  four  regions  of  interest.  Regions 
I,  III,  and  IV  have  a permeability  equal  to  the  permeability  of  free 
space  y which  for  convenience  will  be  labeled  yn . The  problem's 
spherical  symmetry  suggests  that  a spherical  coordinate  system  such  as 
that  shown  in  Figure  3 be  used  in  the  problem  solution. 


6 


Ampere's  law  states 


Vx1T=  T (13) 

and  since  V*B  = 0,  the  induction  B must  be  the  curl  of  some  vector  field  A. 
The  governing  differential  equation  for  A when  homogenous  and  linear  mate- 
rial are  considered  is  from  equation  (11). 


*We  note  that  a distinction  is  drawn  between  the  operator  V Z called  the 
scalar  Laplacian  operator  and  the  vector  Laplacian  operator  designed  by  $ 
The  vector  Poisson's  equation  in  rectangular  coordinates  can  be  treated 
as  three  uncoupled  scalar  equations  as  shown  below. 


0A  = ~e~ 

■v 


- [ A , A,  ,A  l+r[  A , A , A] 
XLA?  yLA  + A-  + AJ  + 


a2A  a2A  32a 

— ~ !~  + — ! 
3x  dyZ  3zZ 


= eJ  + e J + eJ 

xx  y y z z 


where  ^A^  =V  A^  =HJ±  for  i = x,y,z.  However,  if  the  vector  Poisson's 
equation  is  resolved  into  orthogonal  components  in  other  coordinate 
systems  the  differential  operation  mixes  the  components  together  giving 
coupled  equations  as  shown  below  for  spherical  coordinates. 


1 A 2 


. [s_A  9A,  , , A cote 

r L Sr2  r " r2  r r2  Se2  r2  ! 

Ar  2 3A0  2 cot  0 . 2 

2 . 2a  _ , 2 “ 2 ae~  2 e ~ 2 , Q go,  + 

r sin  0 r r r sin  0 J 


A + 

be 


e 2 cot  9 _ 2 _ 

0 2 0 " 2 . . g< 

r r sin  0 


. Sr2 

r dr  ‘ r2 

2 h 2 2 + 
sin  0 r 9 0 

2 d 0 2,2.  ^ 2 

r r sin  0 dip 

dA 

r 

2 cot  0 8AiJ» 

1 PS  + 2 -N 

! . 1 A 

■ ae 

r2  sin  0 *'11 

J 4ar2  r ar 

2 4 2«  A 2 >fi2 

r sin  0 Y r ob 
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3% 


2 2 
r sin  0 


r sin0 


2 cote  dA^ 

~2  — <T 

r sin  0 


= J e + JQeQ  + J . e . 
r r 0 0 ^ i|> 


The  general  expression  in  spherical  coordinates  for  a current  density  is 

■*-rrJr+Ve+V*  (15) 

where  the  X are  the  unit  orthogonal  vectors.  For  stationary  currents 
in  vacuum  the  vector  potential  function  that  satisfies  equation  (14)  is 


where  dv  = elemental  volume  in  the  current-carrying  region 

r'  = distance  between  the  field  point  where  A is  being  determined 
and  dv  at  the  source  point. 

From  equation  (16),  we  see  that  the  elemental  vector  potential  dA 

due  to  a current  element  Jdv  is  in  the  same  direction  as  J.  It  is  well 

known  from  this  that  the  lines  of  the  magnetic  vector  potential  A are 

circles  centered  about  the  coil  or  loop  axis.  The  magnitude  of  X along 

such  a circle  is  constant,  which  means  that  X is  a function  of  the 

spherical  coordinates  r and  0 only.  Therefore,  we  know  in  advance  for  this 

problem  that  A^  is  the  only  component  of  X existing  at  the  field  point. 

The  infinitesimally  thin  band  of  current  shown  in  Figure  2 has  only  an 

azimuthal  or  4*  component,  which  is  a function  of  r and  0,  and  lies  on  the 

boundary  between  regions  III  and  IV  (i.e.  r * R^) . For  this  current, 

equation  (15)  reduces  to 

!,  if  0<0,  or  0>0, 

0 1 2 

(17) 

» if  0 ^ s 0 s ©2 


, if  0jS  0 s 


Therefore,  equation  (14)  has  only  an  azimuthal  component  and  can  be 
expressed  as: 


given  by 

-► 

■ l 

J A 

~ 1 dv 

(16) 

] 

V 

10 


0 (in  regions  I through  IV) 


(18) 


When  the  vector  Laplacian  O is  expanded  in  spherical  coordinates,  equation 
(18)  can  be  written  as 


| 2 [ 1 f cot  6 \ 

3r2  r 9r  r2  9 02  r2  90  r2sin20 


= 0 in  regions  I 

through  IV  (19) 


In  order  to  solve  equation  (19) , it  is  necessary  to  obtain  the  general 

solution  in  regions  I through  IV.  Thus,  by  multiplying  equation  (19) 

2 

by  r we  obtain 


2-2. 
r 9 A 


)k  + + 


92A 


9r 


9r 


..  9A. 

f + cot0  , 

90  9 0 sin  0 


JL 


= 0 


(20) 


Applying  the  method  of  separation  of  variables,  let  us  assume  that  A^ 
can  be  expressed  as  a product  of  two  functions 


A^  = R(r)0  (0)  (21) 

where  R(r)  is  a function  of  r only  and  0(0)  of  0 only.  Substituting 
this  form  of  the  vector  potential  A^  into  equation  (20),  we  have  after 
separation  of  variables 


d2R(r)  2 dR(r)  P(p+l)R(r)  n 

2 ~ 1 “ u 
dr  r dr  r 


(22a) 


d20(0) 

d02 


+ cot  0 


d0(0) , 
d0 


P(P+1) 


1 


sin20  J 


0(0) 


(22b) 


where  the  separation  constant  is  p(p+l)  and  p is  an  integer  from  one  to 
infinity.  The  differential  equation 
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+ cot  0 


has  as  a general  solution 


& 4 

de  L 


p (p+D  - -V  0 
sin  0 J 


0(0)  = 0 (0)  = C Pm(cos0)  + D Qm(cos0)  (24) 

p r p PP 

2 

Comparison  of  equations  (22b)  and  (23)  shows  that  in  equation  (23)  m is 
equal  to  1.  This  requires  that  m always  be  unity.  The  solutions  of 
equations  (22a)  and  22b)  are  then  expressed  as 

R(r)  = R (r)  = A'  rP  + B;  r_(p+1) 

P p p (25) 

0(0)  = 0 (0)  - C P1(cos0)  + D Q1(cos0) 

P p p p p 

The  associated  Legendre  functions  of  the  first  and  second  kind  are  designa- 
ted as  P™(cos0)  and  Q™(cos0),  respectively.  Therefore,  the  general 
solution  of  equation  (19)  in  regions  I through  IV  may  be  formed  from  the 
product  of  the  solutions  in  equation  (25)  which  yields 


A = R(r)0  (0)  = E R (r)0  (0) 
V p=l  P P 


= E/a'  rp+_li\/c  P1 
P=1  \ P r(p+D  / | p P 


(cos0)  + D Q (cos0) 
P P 


In  the  spherical  case,  associated  Legendre  functions  of  the  second  kind 
are  infinite  at  cos0  = +1,  and  thus  cannot  be  included  when  the  region 
under  consideration  includes  the  symmetry  axis.  Therefore,  the  constant 
0^  must  be  set  equal  to  zero,  and  equation  (26)  reduces  to 


where 


?1(vP  + ^+r)  p£(cos0) 


B - B ' C 
P P P 
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BOUNDARY  CONDITIONS 


The  form  of  the  potential  in  each  of  the  regions  (I  through  IV)  is 
determined  from  equation  (27).  These  magnetostatic  vector  potentials  in 
regions  I-IV  are: 

OO 

A = A = V'  fH  v-Pl 


AI  = = E <Aplr  > Pp(cos0) 

p=l 


0 

ah  ~ Aij/n 

00 

- E 

P=1 

A rp+  B-^-l 

P2r  + rCP+l) 

P^CcosO) 

P 

(28) 

i 

II 

r 

Am = \in 

oo 

= E 

P=1 

r b 

A rP  + -P  3 

p3  r(p+l) 

I P^(cos0) 

J P 

II 

in 


Aiv  = ai/>iv  = E 

p=: 


oo 

r p1 

“*  (p+1)  p 
=1*-  J 


(cos0) 


where  for  the  A^  component  B^1  = 0 because  at  r = 0 the  potential  must  be 
finite  and  for  the  A^y  component  A ^ = 0 because  as  r approaches  infinity 
the  potential  must  remain  finite. 

At  each  interface,  the  basic  laws  of  magnetostatics  in  equations  (la) 
and  (lb)  reduce  to  boundary  conditions  on  B and  H that  can  be  used  to 
evaluate  the  six  constants  in  equation  (28).  From  equation  (lb),  the  normal 
component  of  B across  each  boundary  must  be  continuous,  i.e.,  (B^  - "b^) ."n^O 
where  the  quantity  n^  is  the  unit  outward  normal  to  the  surface  . This 
provides  the  following  boundary  conditions  which  must  be  satisfied  by 
the  solution  in  equation  (28)  for  each  region. 


Brl  - Brll  at  r = \ 


Brll  " BrIII  at  r = R2 


BrIII  " BrIV  at  r = R3 
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The  normal  component  of  the  magnetic  field  is  expressed  in  terms  of 
the  vector  potential  as 


Br  = (VxA) r 


B = (sinOA,  ) 

r rsinO  30  ip 


where 


B = VxA  = ~ 2 


r sin0 


dr 


e0r 

3 

33 


e^rsinO 

_3_ 

dip 


r sinO 


(30) 


\ 


However,  since  the  vector  potentials  in  each  region  are  functions  of 
Pp(cos0)  we  can  simplify  equation  (29)  to  constraints  on  A^ 


Aj  = Ajj  at  r = 


Ola) 


'II 


AIII  at  r = R2 


(31b) 


III 


'IV 


at  r = R, 


(31c) 


The  second  set  of  boundary  conditions  is  obtained  from  equation  (12b) 
The  tangential  component  of  H across  each  boundary  must  satisfy  the 
relationship 

x (Hi  - It)  = J (32) 

12  2 1 s 

where  (which  equals  J (0))  is  the  real  surface  current  density  in  the 

limit  of  vanishing  width  between  the  two  regions.  Using  the  relationship 
B = UH,  equation  (32)  may  be  expressed  as 


02 


01 


— - — = J (0) 


1A 


03) 


Referring  to  the  curl  in  equation  (30) , we  can  write  Bg  as 

Be  ■ <?,*>„  - - i [r  Aj;]  (34) 

From  equations  (32),  (33),  and  (34)  the  tangential  components  in 
regions  I-IV  must  satisfy  the  relationships. 


^ 7 (rAII>  + ^ r Sr  <rV  ' ° at  r - R1 


(35a) 


iJ7  r d7  (rAni>  + iJJ  7 ar  (rAn)  = 0 at  r = R2 


(35h) 


^ r ~ (rAIV)  + ^ r hr  (rAIH)  = ^6)  at  r = R3 


(35c) 


The  general  expressions  for  the  potentials  in  each  region  (equation  (28)) 
are  then  substituted  into  the  boundary  conditions  (equations  (31)  and  (35)) 
and  solved  for  the  constants  A ^ and  B^.  There  are  six  algebraic  equations 
with  six  unknowns  and  the  potential  in  each  region  can  then  be  specifically 
determined.  The  six  boundary  value  equations  that  must  be  solved  for  the 
coefficients  are  given  below  (where  the  index  p is  odd  only  and  understood 
to  take  on  values  from  1 to  00  ) . It  is  noted  that  the  current  J^(0)  must  be 


expanded  into  a set  of  associated  Legendre  functions  in  order  to  evaluate  the 

constants  A . and  B . . The  detailed  expansion  is  in  the  next  section 
pl  pl  / \ 

entitled  "Expansion  of  the  Current  ^J^(9)j  1°  Associated  Legendre  Polynomials'  . 


A R P 
pl  1 


■[' 


Ap2RlP  + Vl 


•(P+1) 


] 


(36a) 


[Ap2R2P  + Bp2R2'<P+1>]  ’ [Ap3R2P  + Bp3R2'<P+1)] 


IAp3R3P  + Bp3R3'<P+1>]  ' V[R3"<P+1>] 


(36b) 


(36c) 
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J ■*-  L - 


- t VP+,)  R2  - PBP3S2 


?]4  -J-fV'*1’  "2’ 


- pBo2R2 


-P-2  = 0 


( 36e) 


i~  [ -'vs^]4^  ^p3(P*i)R3p-1  - pbp3'<3'p'2] 


P (cos9) 
P 


The  solution  of  these  equations  to  obtain  B 3 in  terms  of  known  quantities 
is  performed  in  Appendix  A.  In  summary: 


- ~ Jp"  (0)(  [x](p+i)r2(p_1)  -(p)r2  p_2  )+  -J-  Jp'(0)  (p+1)  r2(p  '1) 

i-  1 

. — . — ■ ■ - ■ ■ - ■ — ' ........  ■ ■■-  . — ■ .i  — — , — 

(p1R2~P~2)+  ([’][>]  <P+13R2<p-1))  - ( [z]  Cp3R2‘P_2) 


p3  r/ 1 


where 


-<2p+i>  r /_e_\  ,vil 
[x]-  1 


(■-*-) 


(37b) 


-(P+1) 


[*]- 


([X]R2P  + R2-(P+I)j 


(37c) 


J "(0) 
p 


J »(0)  R P 
P L 


^[X]R2P  + R. 


-(p+1) 


(37d) 


J ' (0)  = — 

P P^cose)*^  '(2p+l) 

p J 


(37e) 


The  numerical  values  for  the  other  five  coefficients  can  be  obtained  from 
the  following  aquations: 


B . - B [Z]  + J " (0) 

p2  p3  L J p 


A - = J ' (0) 
p3  p 


(38a) 


(38b) 


Ap2  = L X1  V 


(38c) 


A = A + B R _(2p+1) 
pi  p2  p2  1 


(38d) 


V ■ Ap3R3<2P+1>  + Bp3 


(38e) 


Since  the  coefficients  A . and  B , can  be  determined  from  equations 

pi  pi 

(37)  and  (38),  equations  (28)  can  now  be  used  to  completely  specify  the 
potentials  A^,  A^,  A^^  , and  AIV  in  regions  I through  IV.  Then  the  normal 
(B^)  and  tangential  (Bg)  components  of  the  magnetic  induction  in  regions 
I through  IV  can  be  determined  by  using  equations  (30)  and  (34),  respec- 
tively. 

In  Appendix  B the  magnetic  vector  potentials  A^  in  the  inner  region 
and  A^jj  in  the  outer  region  are  derived  for  the  infinitesimally  thin 
current  band  in  a homogeneous  medium  of  permeability  y^  (see  Figure  1-B 
in  Appendix  B) . Also,  the  dipole  potential  term  in  the  outer  region 
(r>R^)  for  the  infinitesimally  thin  current  band  is  reduced  in  a special 
case  to  the  dipole  potential  term  for  the  circular  filamentary  current  loop. 
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In  Appendix  C the  coefficients  A ^ (i  = 1,2,3)  and  (i  = 2,3,4) 
for  the  vector  potentials  for  the  present  ferromagnetic  shell  problem  reduce 
to  the  potentials  in  the  two  regions  of  the  simple  current  band  problem 
when  the  permeability  of  the  ferromagnetic  shell  y ^ approaches  that  of 
the  surrounding  medium  y^.  This  shows  that  the  solutions  of  the  above 
ferromagnetic  current  problem  have  the  correct  mathematical  form. 

EXPANSION  OF  THE  CURRENT  (j^(0))lN  ASSOCIATED  LEGENDRE  POLYNOMIALS 

Any  function  that  can  be  expanded  using  a Fourier's  series  in  a given 
interval -l<yf<l  can  also  be  expanded  in  a series  of  associated  Legendre 
polynomials  in  the  same  interval  using  similar  methods.  The  associated 
Legendre  functions  are  defined  by  the  equation 


m 


(y1)  = (1  -y  ) 


.2^  £VP') 


dy 


i m 


(39) 


where  y*  = cos0  and  -liy'sl.  Also,  the  function  P^(y'),  which  is  valid 
whatever  the  range  of  the  variable  y1,  is  defined  as 


P (y1)  = 
P 


2pp! 


dy 


IP 


(y,2-i)p 


Let  us  assume  the  expansion  is  similar  to  Purczynski ' s^ 


•ye>  = £ J Kp  Pp(cos0) 
p=i 


(40) 


(41) 


J 

1 

3 


■ J 


The  coefficients  K are  determined  from  equation  (41)  by  multipling 

1 P 

both  sides  by  P (cos0)  and  integrating  over  sin0d0  from  0 to  n . 

P 


K 


122+1). 

2p (p+1 


>*/ 


J|/,(0)  P Ccos0)  sin0d0 
r0  V p 


(42) 


[1 

fl 
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The  identity 


/ 

'*'0 


p",'  (co3e)  p;<'°*e)  s i„ede  - Vp 


T « 


(43) 

was  used  to  determine  K (equation  42).  Setting  J, (0)  to  a constant  J 

P V 

we  can  express  equation  (42)  as: 

•7T 


K =2£±1 

p 2p (p+1) 


/•TT 

I ' 

^ r\ 


( cosQ ) sinddd 


(44) 


Since  the  current  in  this  problem  extends  from|^ — aj to  + a j , 
see  Figure  4,  the  expression  for  K (equation  44)  may  be  written  as 


K = 2P+1 

p 2p(p+l) 


r < 

- a 


(cos0)  sin0d0  + 


/"**  p1 
A p 


(cos0)sin0d0 

(45) 


By  noting  that 


(-cos0)  * (-1)^  ^ P^(cos0) 
P P 


(46) 


it  follows  that  the  associated  Legendre  functions  Pp(cos0)  are  even  func- 
tions with  respect  to  cos0  when  p is  odd.  The  expression  for  K may  be 
simplified  to 


K 


2p+l 


I, 


P “ p(p+l)  J it  Pp(cos0)  sin0d0.  (47) 

2 " “ 


when  p is  odd  (K^  = 0 when  p is  even)  after  utilizing  the  symmetry  of  J(0) 
in  Figure  5.  By  changing  the  variable  in  equation  (47)  to  y 1 = cos0,  the 
integral  for  may  be  written  as 


K = 2P+1 

p p(p+l) 


>.  sina 

/ Pj(u')dy' 

•'0 


(48) 
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The  values  of  the 
9 

Purczynski  and  are: 


coefficients  K 

P 


for  p odd  were  worked  out  by 


3 

Kj  = (sinacosa  + a) 

(49a) 

7 

64 

sinacosa(10  sin2a-9)  + aj 

(49b) 

11 

256 

£sinacosa(56  sin^a-70 sin2a+15)  + aj 

(49c) 

* 


^ ^sinacos«(13728  sin^a-23408  sin^a 


7 32768 


+ 11060  sin  a-1330)  + 50 


(49d) 


K9  = 327680  £ slnacoaa( 311168  sin®a-679536  sin^a 
+ 488488  sin4a-128590  sin2a+8715)  + 245  a] 


(49e) 


J 

} 


I 


i 


r 


Figure  5 - Coil  Current  as  a Function 
of  the  Polar  Angle  Theta 

FERROMAGNETIC  SPHERE  SURROUNDED 
BY  A COIL  OF  FINITE  WIDTH 

GENERAL  SOLUTION 

We  now  proceed  to  solve  the  boundary  valve  problem  of  a solid  ferro- 
magnetic sphere  of  radius  R^  and  homogeneous  permeability  y^  surrounded 
by  a current  band  of  finite  width  having  inner  radius  R£  and  outer  radius 
R_  as  shown  in  Figure  6.  A constant  current  density  is  assumed.  A linear 
relationship  between  B and  H is  assumed.  Regions  II,  III,  and  IV  have 
homogeneous  free  space  permeability  designated  as  y The  permeability 
of  the  conducting  coil  in  region  III  is  assumed,  also,  to  be  equal  to 
The  geometry  of  the  problem  suggests  spherical  symmetry  as  in  the  previous 


problem. 

The  partial  differential  equation  that  governs  this  problem  is  again 
equation  (11), 


where 


$T=-yT 

VVr,0) 


(50) 
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CONDUCTOR 
CURRENT  DENSITY 

7-9.1.  J.i.  (r,  o) 


FERROMAGNETIC  MATERIAL 


Figure  6 - Ferromagnetic  Sphere  Surrounded 
by  a Coil  of  Finite  Width  (yz  Plane) 


As  before,  it  is  necessary  to  solve  equation  (50)  and  use  the  appropriate 
boundary  conditions  to  evaluate  the  constants.  Equation  (50)  reduces  to 


where 


O A,„  = - *Vr*9> 


J,  (r,0)  = J,  (0)  in  region  III  if  and  R^rsR^ 


J^(0)  = 0 if  0<0^  or.  0>©2  for  all  r(0^  and  ©2  are  defined  in 


Figure  6) . 


Expanding  the  vector  Laplacian  in  equation  (51)  in  spherical  coordinates 
results  in  the  expression: 


- 


k_ 


+ i + I_ + ^ ajl_ 

Sr^  r &r  r2  902  r2  90  r2sin20 

(52) 

- pJ^(0)  in  region  III 
* 0 in  regions  I,  II,  and  IV 

Thus,  the  solutions  in  regions  I,  II,  and  IV  are  solutions  of  Laplace's 
equation  and  are  obtained  by  the  method  of  separation  of  variables  as  in 
the  previous  problem.  The  potential  solution  in  region  III,  however,  is 
the  solution  of  Poisson's  equation. 

The  potentials  in  regions  I,  II,  III,  and  IV  are: 

ai  ■ v <vp>  pJ(c°s6)  ™ 


AII  ■ Vi  ■ £ ( VP  + T^i)  ) pJ<“sS>  (53W 

p=l  ' r ' 

AIII  = \lll  = £ (Ap3rP  + r(p+l)  ) Pp(cos6)  + 


/ Vr  Kn  \ 1 

p«\(p-2)(p^3)  ) V“se) 

(53c) 

P = 1 

AIV  Vv‘  (r(p+i)  )pp(cose) 

(53d) 

Equation  (53c),  the  azimuthal  component  of  Poisson's  equation  in 
spherical  coordinates,  is  solved  in  detail  in  Appendix  D.  When  the  current 
is  expanded  in  associated  Legendre  functions  we  have 


j 


* t 
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J i (r , 0 ) = J(0)  = J £ K p1(cosG)  (54) 

* p-1,  P P 

where  only  the  odd  coefficients  contribute  since  the  even  coefficients 
are  zero.  The  coefficients  for  the  current  were  derived  previously  (see 
equations  (49a  through  49e)). 

BOUNDARY  CONDITIONS 

This  solution  must  satisfy  the  boundary  conditions  developed  in 
equations  (29)  through  (33).  Because  there  are  no  surface  currents  on 
the  boundary  between  the  regions  in  this  ferromagnetic  problem,  Jg  is  zero 
and  the  tangential  boundary  conditions  for  regions  I-IV  become 


— — (rATT(r,0))  = — ^ (rAT(r,0))  at  r = R. 

D2  r dr  II  ’ y^  r Sr  I 1 

(55a) 


i £ <rAin(r,9)) 


h (rAn(r,6))  .t  t - «2 


(55b) 


y.  r Sr  (rAlv(r,0))  “ y r Sr  'rAIII(r,0))  at  r ' R3 


(55c) 

The  general  expressions  for  the  potentials  in  each  region  (equation 
(53))  are  then  substituted  into  the  boundary  conditions  (equations  (31)  and 
and  (55))  and  solved  for  the  Apiis  and  Bpi's*  T^ere  are  s*x  algabraic 
equations  with  six  unknowns  and  the  potential  in  each  region  can  then  be 
specifically  determined.  The  six  boundary  conditions  that  must  be 
solved  for  the  coefficients  are  (where  the  index  p is  odd  only  and 
understood  to  take  the  values  1 to  «) : 


SSI 


A -R/ 
pi  1 


A „R  P + B „R  "(p+1) 


(56a) 


f 


A.,R,P  + B_„R  ~(p+1)  = A.,R,P  + B_  _R_-(p+1)  + 


p2  2 p2  2 


p3  2 p3  2 


(p-2)(p+3) 


(56b) 


A „R  P + B R -(p+1)  + .-2-  y,-^r.  = B ,R  -(p+1) 
p3  3 p3  3 (p-2) (p+3)  p4  3 


(56c) 


i-  Ap2(P+1)R1p-1  -(p)Ep2<t1-t’-2  j - ^ [ap1(P+1)»1'’-1  j 


(56d) 


i ^p3<p+i)e/-1  -(p)Bp3R2-p-2  4- 


T.  [ VP+1)R/'1  -(p)Bp2R2"p"2] 


(56e) 


57  [-(p>Bp4R3'P'2]  ■ 17  [AP3<P+1)R3P'1 

3P2JR3ICp  1 
(p-2) (p+3) J 


•<p)Bp3R3'P'2  + 


(56f ) 


The  mathematical  solution  for  B^  in  terms  of  known  quantities 
obtained  in  Appendix  E and  is  given  by 


( , 3y0JR0K 

Bp3  ‘ -Kp'  <^»R2<P'1)  * 7P2K^3) 


(57a) 


+ [x]  (p+DR/'V*  - (p)R2"p"2Kp' ' 


j (-p)R2'P^2  - [Z][x]  (p+l)R2P_1  +[z]  (p)r2"p"2J 
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~(2p+l)  i+/  JL 


(-S) 


U„JK  R_  * 

2 P 3 

(p-2) (2p+l) 


(57c) 


-(P+1) 


[Z]  '("[x]  R/ + r2-(p+1)  )' 


(57d) 


y0JR0  K 

2 -2  P + K ' R P 
'*  = (p-2) (p+3)  + P K2 

P |[x]  D P I d "(P+D  ' 


R2  + R2 


(57e) 


The  numerical  values  for  the  other  coefficients  can  be  obtained  from  the 


equations 


aP3  = y 


(58a) 


Bp2  " Bp3  0]  +KP” 


(58b) 


" W K 


(58c) 


A = A + B R "(2p+1) 
pi  p2  p2  1 


(58d) 


,,  ...  U_JR.(p+3)K 

B = A _R_ ^2p  + B + 2 .1.  ..'F 

pA  p3  3 p3  (p-2) (p+3) 


(58e) 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  method  of  separation  of  variables  has  been  applied  to  determine 
the  magnetic  field  of  the  systems  consisting  of  an  infinitesimally  thin 
spherical  coil  around  (outside)  a ferromagnetic  spherical  shell  and  a 
spherical  current  band  of  finite  width  around  a solid  ferromagnetic  sphere. 
The  resulting  formulae  can  be  used  in  the  analysis  of  magnetic  induction 
of  ferromagnetic  bodies  due  to  current-carrying  coils. 

The  magnetic  vector  potentials  in  the  inner  and  outer  regions  are 
derived  for  the  infinitesimally  thin  current  band  in  a medium  of  homo- 
geneous permeability.  The  dipole  term  for  the  potential  on  the  outer 
region  for  the  infinitesimally  thin  current  band  is  reduced  in  a special 
case  to  the  dipole  potential  term  for  the  circular  filamentary  current 
loop.  It  is  also  shown  that  the  vector  potential  for  the  ferromagnetic 
shell  surrounded  by  a infinitesimally  thin  current  band  reduce  to  the 
potential  in  the  two  regions  of  the  simple  current  band  problem  when 
the  permeability  of  the  ferromagnetic  shell  approaches  that  of  the 
surrounding  medium. 

It  is  recommended  that  the  solutions  derived  in  this  report  be 
programmed  on  a digital  computer.  The  resulting  calculations  should  then 
be  compared  to  solutions  obtained  by  various  numerical  methods  in  order 
to  validate  the  numerical  methods  and  calculations.  There  are  plans  to 
implement  these  recommendations  during  the  fiscal  years  1979  and  1980. 
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APPENDIX  A 

CALCULATION  OF  COEFFICIENTS  OF  THE  VECTOR 
POTENTIALS  FOR  A THIN  COIL  SURROUNDING 

A FERROMAGNETIC  SPHERICAL  SHELL 

In  this  appendix  the  coefficients  are  derived  for  the  vector  poten- 
tials in  regions  I-IV  for  a ferromagnetic  spherical  shell  surrounded  by  an 
infinitely  thin  current  band.  For  a detailed  discussion  of  the  ferro- 
magnetic problem  see  the  section  in  the  text  of  the  report  entitled  "Thin 
Coil  Surrounding  a Ferromagnetic  Spherical  Shell".  The  magnetic  vector 
potentials  in  each  region  are  given  by: 

OO 

\l  " ^ ^AplrpjPp(cos0)  (A-la) 

OO  / B \ 

A»n ' P?i  ( vp + -feu) p («•«  «-ib> 

A<I>III  " ^ ^Ap3rP  + r(p+l)  ) pp(cos9)  (A-lc) 

Aipiv  ■ z pJ<cos6)  <A-ld) 

The  coefficients  (Ap^  and  Bp^)  in  equations  A-la  to  A-ld  are 
obtained  by  substituting  these  equations  into  the  boundary  conditions 
(equations  (A-2a  to  A-2f)) 


i 


(A-2d) 


1_1.S 

y2  r dr 


(rAu>  + 57  7 »7  (tV  ‘ 


0 at  r 


- ^ 7 »\  «AIII>  + i2  'r  '*  <rtII>  * 0 at  r ■ R2  (A-2e) 

' "r  3r  <rflIV>  + r hr  <rAIII>  ' J<e>  at  r ’ R3  <A'2f) 


After  appropriate  substitutions  of  equations.  A-la  to  A-ld  into  equations 
A-2a  to  A-2f,  the  following  boundary  value  equations  are  obtained. 


(ApiV> 


I 

[ 


A „RJ 

p2  1 


E? _ \ 

f i \ ! 


(p+1) 


(A- 3a) 


(A- 3b) 


(A-3c) 


- 57  - <P>  VRl"”'2] 

+ J-  ^Api<p+1)RIP'1  J -0  (A-3d) 

- ^ [Ap3<p+1)R2'>'1  * (p>  Bp3R2"P'2] 

+ ^ [ Ap2  <p+1>r2P'1  - <P)  Bp2r2"P“2  ] ' 0 (A-3e) 
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- - I (-P)Bp4R3 


-P-2 


] +k  [ 


Ap3<p+l)R3(p-1>  -<P>Bp3R3-»-2 


V6> 


P (cos0) 
p 


(A-3f) 


These  algebraic  equations  provide  six  simultaneous  equations  with 

six  unknowns,  and  they  can  be  solved  for  the  coefficients  A . and  B , by 

Pi  Pi 

tedious  algebraic  manipulation. 

Solving  equations  (A-3a  and  A-3d)  algebraically  results  in  the 

solution  for  A _ in  terms  of  B , 
p2  p2 . 


V [x]  V 


(A-4) 


where 


[X] 


TJ  -(2p+l) 
R1 

[>♦(*) 

1 9 

1 

y2y 

1 

The  solution  for  Ap^  from  equation  (A- 3a)  is 


Apl  - V + Bp2Rl 


-(2p+l) 


(A-5) 


Also,  solving  equations  (A-3c)  and  (A-3f)  algebraically  for  B^  gives  the 
expression  for  Ap^  it*  terms  of  known  quanitities  after  equating  the 
functions  for  B^  to  B^  from  each  equation.  Thus, 


A = J'(0) 
p3  p 


(A-6) 


where 


Jp'(0) 


^lJp(9> 


P 1(cos0)  R,(p_1)(2p+1) 
P 3 
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The  algebraic  solution  for  B ^ from  equation  (A-3c)  is 

B , = A ,R  (2p+1)  + B 
p4  p3  3 p3 


(A- 7) 


The  following  expression  is  obtained  for  B^  from  equation  (A-3b) 

after  substituting  the  expressions  for  A _ (equation  A-6) , and  A „ (equa- 

P -J  p2 

tion  A-4) . 


Bp2  - Bp3  M + Jp<e> 


(A-8) 


where 


[*> 


R„ 


•(p+1) 


([X]  R. 


,p  + v(p+1)) 


and 


J'(9)R?P 
ji.(e)  E L 


~(M 


R P + R„ 


(p+Dj 


p3 


R2'P'2V^tZKx]<P+l)»2<P'1)Vr/'tz]  pE 


)"M‘ 


-p-2 


7 


(A-9a) 


where 


[X] 


- R, 


(2p+l) 


H) 


(A-9b) 


i 

J 

I 


The  mathematical  solution  for  B^  in  terms  of  known  quantities  is  obtained 

from  equation  (A-3e)  by  substituting  the  previously  obtained  expressions  ; 

for  A « (equation  A-6) , A 0 (equation  A-4) , and  B 0 (equation  A-8) . 

P-3  pz  pz 

- TT^  Jp*(S>(f  x3  (p-|-»k2(p~1)  - (p)r2~p~2  j + jL  Jp'(e>  (p-h)r2(',~')  | 


Jl 


0 

II 

0 
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r2~(p+1) 

[ x]  R2P  + R2"(p+1)' 


V,(0) 


j '(e)  r p 

P 2 


^[x]  R2P  + R, 


-(p+l)> 


(A-9c) 


(A-9d) 


H J (0) 

J ' (0)  = 7 ^ 7 T7T 

P Pi(cos0)R  P"i;(2p+1) 


(A-9e) 


After  the  numerical  value  for  is  calculated  on  the  computer  for  a 
specific  problem,  the  numerical  values  for  the  other  coefficients  can  be 
obtained  from  the  following  equations: 


6^2  = Bp3  [ z]  + Jp"(9),  (see  equation  A-8) 


A^  = (0),  (see  equation  A-6) 


A^  = [x]  B^2,  (see  equation  A-4) 


= A^  + Bp2R1  ^P+l),  (see  equation  A-5) 


(A-lOa) 


(A- 10b) 


(A- 10c) 


(A-lOd) 


Bp4  * Ap3R3^P+^^  + Bp3  ^see  e9uat*on  A-7) 


(A-lOe) 
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APPENDIX  B 


DETERMINATION  OF  POTENTIALS  FOR  INFINITESIMALLY 
THIN  CURRENT  BAND  AND  INFINITESIMALLY  THIN  COIL 


0 

0 


In  this  appendix  the  potentials  A^^  in  the  inner  region  and  A^jj  in 
the  outer  region  are  derived  for  the  infinitesimally  thin  current  band  in  a 
homogeneous  medium  of  permeability  (see  Figure  1-B) . Also,  the  potential 
in  the  outer  region  (r>R^)  for  the  infinitesimally  thin  coil  is  shown  to 
reduce  to  the  potential  of  the  circular  filamentary  current  loop. 

The  potentials  in  the  inner  region  A^  and  outer  region  of  the 

infinitesimally  thin  current  band  prbblem  are  solutions  to  the  vector 
Laplace's  equation  $A  = 0.  These  solutions  can  be  expressed  as: 


V = S > Vcos9) 

p=l 


4 CO  / B \ 

= £ ( ) P1(cos0) 

p=l  \ r(p+D  } P 


(B-la) 


(B-lb) 


The  coefficients  A^  and  B^  are  determined  from  the  boundary  conditions  of 
the  problem.  After  algebraic  manipulation  such  as  with  equations  (31)  and 
(35)  in  the  text,  the  boundary  conditions  for  the  normal  component  of  B 
and  the  tangential  component  of  H become: 


A\pl  AipII  at 


r = R, 


r dr  (rAII)  r dr  (rV 


(B-la) 


J^(6)  at  r = R1 
(B-2b) 


Substituting  the  expressions  for  A^  and  A^II(equations  B-la  and  B-lb)  into 
the  boundary  value  equations  (equations  B-2a  and  B-2b)  provides  us  with  the 
following  algebraic  equations  for  the  coefficients: 


A R P = B ?R  (p+1) 
pi  1 p2  1 


(B-3a) 


u; -isfry  > •;*. 


- . .. 

...  ....  • t&SSSeiffixiSjMk*. T 


i-[  (-P)Bp2VP'2] 


(B-3b) 


+ [APi<p+i>Rr1]  ■ jp<e)/pp<cos°) 

where  3^(0)  is  the  p1"*1  term  of  J^(0)  (see  equation  41).  These  equations 
are  solved  for  A^  and  B^  hy  simple  algebraic  manipulation. 


Apl  Bd2R1 


-2p-l 


(B-4a) 


1^(0)  / P_(cos0) 
p2  R1~P'2(2p+l) 


’Ap1 


(B-4b) 


The  potential  A^^.  is  determined  by  substituting  the  expression  for 

B _ (equation  B-4b)  into  equation  (B-lb) . 

P^ 


>11 


rVp(e>/pi 


p=i 


r (p+l) (2p+i)r 


(cos0)  -1  n 


(cos0) 


(B-5a) 


^ lJiJKnpp<“s0)Ei 
p=l  (2p+l) 


(B-5b) 


where 


and 


J (0)  = JK  p1 (cos0) 
P P P 

jkp  ■ jf&rlo 


(B-5c) 


(B-5d) 


A^j^  is  the  potential  in  region  2 which  is  external  to  the  coil  (see 
Figure  1-B) . 
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Now  the  solution  for  (equation  B-5b)  can  be  reduced  to  that  for 

the  filamentary  coil  and  the  answer  compared  with  that  for  the  same  problem 
from  a standard  text  (Jackson, ^ page  144,  equation  (5.46)). 

The  coefficient  JK^  for  the  expansion  of  the  current 


J .(0)  JK  p]- (cose) 

r p P P 


(B-6) 

for  the  current  filament  is  obtained  from  the  equation  (see  equation  B-5d) : 


■ i S ’ 4 

■S  :>  4* 


II 

D 

fl 


JK 


2p+l 

2p(p+l) 


X 


COS(TT)  T i 

— 6 (cos  9) P (cos0)d(cos0) 
Ri  P 

cos(O)  1 


(B-7) 


_ 1_  / 22±1_  \ pl 

R1  \ 2p (p+1)  / P 


(0) 
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where 


- (2p+l) 

2p(p+l) 

p 

R-  p 

1 

J(Q)  = 6(cos0) 

R1 


The  special  value  of  the  associated  Legendre  function  is  calculated 
from  the  expression  (where  m = 1) 


, (p+m,  even) 
» (p+m,  odd) 


(B-8) 


This  expression  may  be  derived  from  the  series  expansion  for  the 
Legendre  function  (see  Smythe) ^ and  letting  cos0  = p1  approach  0, 

Substituting  the  expression  for  JKp  (equation  B-7) , into  equation 

(B-5b),  results  in  the  expression  for  the  vector  potentials  A,„  in  the 

(JjII 

external  region  for  the  filamentary  coil. 


\ll  = 


y,lP1(cos0) 
2 

(2p) (p+1) 


P+1 


pj<0) 


(B-9a) 


Now  equation  (B-9a)  must  be  multiplied  by  ~ to  change  from  SI  units 
to  Gaussian  units,  and  p^  must  be  set  equal  tp  1 as  p^  is  the  free  space 
permeability  which  equals  1 in  the  Gaussian  system. 

After  making  these  substitutions  we  have 


4,  £ i\  f1  ,1(0) 

~ \ 1 I f 


(B-9b) 


Now  taking  the  first  term  (p=l)  for  use  in  the  comparison  with  the  result 
in  Jackson^  we  have 
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P^(cos0) 


(B-9c) 


This  expression  agrees  exactly  with  the  expression  in  Jackson10  for  the 
first  term  in  the  series  where  R = a in  his  notation.10  The  magnetic 
fields  far  from  the  loop  for  the  p=l  term  are  dipole  in  character. 


term  for  the  current  band  problem  also  reduces  exactly  to 


Jackson’s  term  with  n=0 


(Jackson, 


10 


page  144,  equation  5.46)  for  the 


potential  inside  of  the  loop. 


APPENDIX  C 


REDUCTION  OF  THE  MAGNETIC  VECTOR  POTENTIAL 
FOR  A THIN  COIL  SURROUNDING  A FERROMAGNETIC 
SPHERICAL  SHELL  TO  THAT  OF  A THIN  COIL 
IN  FREE  SPACE  WHEN  IN  THE  LIMIT  y 
EQUALS  u2  1 


In  this  appendix,  the  coefficients  Ap^>  Ap2’  Ap3’  Bp2’  Bp3’  and  Bp4 
for  the  potentials  are  evaluated  for  the  system  consisting  of  a ferro- 
magnetic shell  with  permeability  surrounded  by  an  infinitesimally  thin 
current  band  in  a homogeneous  medium  with  permeability  y^  in  the  limit  as 
^1  = y2*  T^ese  coefficients  are  utilized  in  equation  (28)  in  the  section 
of  the  report  entitled  "Thin  Coil  Surrounding  a Ferromagnetic  Spherical 
Shell".  The  variables  are  defined  in  Figure  2 located  in  the  text  of 
this  report.  When  y^  is  set  equal  to  y2  the  problem  reduces  to  that  of 
finding  the  potentials  in  the  two  regions  of  a simple  current  band  (see 
Figure  1-B  in  Appendix  B) , since  the  ferromagnetic  shell  will  now  have 
a permeability  y^  equal  to  that  of  the  homogeneous  medium  with  permeability 

yr 

In  this  limit  the  coefficients  should  assume  the  following  form: 


>1  P2  P3 


(C-la) 


p2  P3 


(C-lb) 


and  where  A^  and  B ^ should  reduce  to  the  coefficients  for  the  potentials 
in  the  two  regions  for  the  spherical  band  problem  (see  Appendix  B).  If 
the  coefficients  assume  this  mathematical  form  it  will  prove  that  the 
mathematical  form  of  the  coefficients  for  the  spherical  shell  surrounded 
by  a thin  current  loop  are  mathematically  correct. 


The  mathematical  solution  for  B^  in  terms  of  known  quantities  was 
derived  in  Appendix  A and  was  reported  in  the  text  of  this  report  (see 
equation  (37a))* 


The  expression  for  J equals  J^(0)  (see  equation  (38b)  in 

the  text  of  this  report).  The  expression  for  (see  equation  (38a)  in 
text  of  this  report)  is: 


V ■ V[z]  + Jp(e> 


The  expression  for  when  can  be  expressed  as: 


(C-7) 


(Uj-Mj) 


= limit 


Ik,  ll^ 


/ Jn(0)  R9P  V 

limit  / 2 7— , . . \ 

1. . ' [*]  *2P  + V<pH7 


(C-8) 


(yi=y2> 


where 

limit 

0]  - 

| li 

[ xj-*-  00 

[*]■ 

and 

limit  J 

< j_'(9)  pJ 

r±-  — L 

n 

w—  \ 

([x]  R2p  h 

1 *y 

and 

| li 

Bp3  J 

<YV 

- 0 

-(p+1) 


([X]  R. 


2P  + r2-<p+i>) 


The  expression  for  A^  is  (see  equation  (38d)  in  the  text) 


Apl  " Ap2  + Bp2Rl 


(2p+l) 


(C-9) 


t 


The  expression  for  A when  p,  = p„  is 

pi  11 


Pi 


II 

N> 

+ 

V 

=P2)  (Pj* 

y2) 

(p^l 

J2> 

- (2p+l) 


Pi 


- j;«) 


W 


where  B 


p2 


= 0 


(y-^Pj) 


The  expression  for  (see  equation  (38e)  in  text)  is: 


(C-10) 


B = A _R  (2p+1)  + B 
p4  p3  3 p3 


(C-ll) 


The  expression  for  B ^ when  p^  = p2  is 


Bp4 

s 

Ap3 

R (2P+1)  + g 

3 p3 

(Mr 

y2) 

- <v 

p2) 

(pr 

(C-12) 


*P3  | 

R (2P+D 

R3 

V | ' 

. W- 

„ (2p+l) 

R3 


where  B „ | * 0 . 

P3 


(yj-yj) 
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This  means  that  in  the  four  regions,  the  potentials  used  in  equation 


(28)  of  the  report  are 

or 

V = £ (AplrP)  PJ(cos6)  (C-13a) 

P-1 

00  / B \ 

V'E  */  pJ(cose)  (c-13b) 

P=P  r ' 

OO  g 

\m-  E ( V"  + Pp1<'‘>se>  (C-13c) 

00 

A^IV  £ | r(p+l) j pp(cos0)  (C-13d) 


(C-14a) 


(C-14b) 


These  are  the  solutions  for  the  potentials  of  the  current  band  in  a region 

of  space  with  homogeneous  permeability  (see  equations  (B-la  and  B-lb). 

We  now  have  the  solutions  for  the  two  potentials  in  regions  I and  II 

(A^i  jj  and  A^IV>  respectively)  for  the  simple  current  band  problem. 

This  indicates  that  the  form  of  the  coefficients  A . and  B . are  mathemati  - 

pi  pi 

cally  correct. 
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The  mathematical  expressions  for  A , and  B . will  be  evaluated  in 

pi  p4 

the  limit  as  These  values  will  then  be  compared  with  the  coeffi- 

cients A^  and  Bp2*  respectively,  for  the  two  regions  of  the  current  band 
problem  (see  Appendix  B) . B ^ (see  equation  (38e)  in  text)  is: 


B , = A R_(2p+1)  + B 
p4  p3  3 p3 


When  \*i  = U2>  Bp4  is: 


= [AP3  I |R3(2P+1)+Bp3 


(p1=m2) 


V • *‘lJp(e)J 


'Pl(cos0)R  "p  2(2p+l) 

P J 


‘W 


The  form  of  A^  I is: 


(VU2> 


(M1^2> 


<W 


-2p-l 


(C-15) 


(C-16) 


(C-17) 


(y1=y2) 

[ ^i=^2>J 

(Wl-U2) 

where 

AP3 

= J'(6) 

P 

r 

(P1! 

CN 

ll 

r 

and 

BP3 

1 ‘° 

The  mathematical  expressions  for  A , and  B . (equations  C-17  and  C-16, 
respectively)  for  the  ferromagnetic  spherical  shell  surrounded  by  a thin 
current  band  in  the  limit  as  p^  = are  the  same  as  the  coefficients 
A^  and  (see  equations  B-4a  and  B-4b) , respectively,  for  the  vector 
potentials  in  the  regions  of  the  current  band  in  free  space  (see  Appendix 
B) . It  is  noted  that  when  making  the  comparison  must  be  set  equal  to 
For  comparison,  the  coefficients  for  the  current  band  problem  are: 


A n = B „R 
pi  p2  1 


ylJp(6)  /vcos9> 

R1"p‘2(2p+1) 


(C-18a) 


(C-18b) 


and  the  coefficients  for  the  ferromagnetic  shell  problem  with  p^  = \i^  are: 


-2p-l 


(C-19a) 


yiJp(9)  / vc°89> 

R3“p_2(2p+1) 


(C-19b) 
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appendix  d 


SOLUTION  OF  THE  AZIMUTHAL  COMPONENT  OF  THE  VECTOR 
POISSON'S  EQUATION  $ A^  = IN 

SPHERICAL  COORDINATES 


For  the  benefit  of  the  reader,  the  azimuthal  component  (ip  component) 

of  the  vector  Poisson's  equation  in  spherical  coordinates  as  originally 
> 9 

derived  by  Purczynski  is  presented  in  detail  as  related  to  the  problems 
addressed  in  this  report.  Because  of  the  spherical  symmetry  of  this 
problem,  only  the  ip  component  of  the  vector  Poisson's  equations  is  needed. 
Following  in  an  outline  of  Purczynski ' s development,  the  general  form  of 
the  component  of  the  current  is : 


J,(r,9)  = Jrq~2  ^ K P1(cos0) 
V p=i  P P 


(D-l) 


where  the  case  q = 2 is  of  primary  concern  in  this  work.  The  general 
component  of  Poisson's  equation  is  written  as 


* \ , 2 


r 9 r 


2 2.  2 .q2 

r sin  0 r 90 


i 32A.  q 9A. 

. JL  1 + cot9  \p 

2 *02  2 
r 00  r 90 


"^2JiJ)(r’0)  *n  current  region 


otherwise 


(D-2) 


Multiplying^^  by  r'  and  substituting  the  general  expression  for 
J, (r,0)  in  equation  (D-l)  gives  the  expression 


+ cotO 


_A_+ 

9,2  »ln26  392 


- y0Jr^  Z K p1(co®0) 


which  is  assumed  to  have  a solution  of  the  form 


(D— 3) 


A - A - £ f(r)P1(cos0) 

P 

P“1 

By  separation  of  variables,  equation  (D-3)  becomes: 


d2  fp^(cosO)l  d |P^(cos0) 


cot  0 + 


(D-4) 


p(p+l)  - 


— — y~~\  P1 

sin  0 J P 


(cosO)  = 0 


(D-5a) 


r2  + 2r  - P(P+Df  = ~ P2  JrqK 

dr2  dr  2 P 


(D-5b) 


where  the  separation  constant  p(p+l)  is  an  integer.  Differential  equations 
of  the  form 


+ p^x*  lx  + = 

dx 


(D-6) 


have  no  general  solution.  The  homogeneous  part  of  equation  (D-5b)  has  the 


2 d f , . df  , ....  . 

r — 2 + 2r  dr  " P(P+Df  * 0 
dr 


(D-7) 


Equation  (D-7)  can  be  simplified  by  substituting  r ■ e* . Making  the  follow- 
ing substitutions  in  equation  (D-7) , 


df  _ df  dt  _ df  / 1\ 
dr  " dt  dr  * dt  \ r / 


(D-8) 


d2f  d2f  /dt\?  df  djt_ 
dr2  dt2W  dt  'drS 


r 


f 


1 

d"f 

-£l/k 

\ +df 

/ i_  \ 

1 

dr2 

l J + dt 

( rO 

where 


— = - , 4»-^,  (Note:  i n(r)  = t) 

dr  r dr  r 


yields 


-P(p+l)f  = 0 

dt^  ac 


(D-9) 


This  has  the  well-known  mathematical  solution: 


f - A ept  + B e"(p+1)t 
1 p p 


(D-10) 


which  after  substitution  of  r = e has  the  form 


f.  = A rH  + 
1 P 


B 

E_ 

.(P+1) 


(D-ll) 


Thus  the  general  mathematical  solution  to  the  homogeneous  part  of  the 
azimuthal  component  of  Poisson's  equation  is 


A1  = A* 


00  . R v 

E Arp+  -Ay-V1 
p-1  \ p r(p+1)  / p 


(cos9) 


(D-12) 


which  is  the  solution  to  the  azimuthal  component  of  the  vector  Laplace's 
equation. 

The  solution  to  the  inhomogeneous  equation  (T)-5b)  will  now  be  investi- 
gated. It  is  assumed  that  the  mathematical  solution  to  equation  (D-5b) 
has  the  general  form 

f2  » dpeqt  (D-13) 

for  p i q.  By  algebraic  manipulation  after  substituting  into  equation 
(D-5b)  we  have: 


a . 


p (p-q) (p+q+1) 


. (p  + q) 


(D-14) 


If  we  assume  that  the  mathematical  solution  to  equation  (D-5b)  has  the 
general  form  for  p = q 


f2  ' 


(D-15) 


It  is  found  after  substituting  f ^ into  equation  (D-5b)  that 


- U,  J K 

d - 9 , (p-q)  (D-16) 

9 2q+l 


Since  A equals  f (r)P^ (cos0) , the  general  mathematical  solution  to 


P-1 


the  inhomogeneous  equation  (D-5b)  has  the  form 

K_P^(cos0)  K_P^(cos0)S,n(r)^ 


L p^q 


2q+l 


The  total  general  solution  to  equation  {D-2)  consists  of  the  sum  of  the 
homogeneous  and  inhomogeneous  equation. 


A . »'  ♦ A"-  *;  + V ■ E (a/  + V_,)  P J 


(cos0) 

(D-18) 


+ PoJ r ■ 


K p1  ( cos0 ) K P1(cos0)«,n(r)| 

J3_a 


T P-P 

p^c  (p-q)  (p+q+D 


P*q 

P-1 


2q+l 


For  the  problem  worked  in  this  report  q - 2 and  Kp  - 0 for  even  p.  Thus, 
Kj  is  0 and  the  term  involving 

2 


-M2Jr 


[■ 


K.P^  (cos6)  fcn(r)"l 

,j 
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in  equation  (D-18)  does  not  contribute.  In  this  case  (D-18)  reduces 

to  (see  equation  (53c)  in  text  of  report  where  A * A „ and  B = B _)  • 

p p3  P p3 


A = A = y 
III  (JiIII 

P-1 


V 


rp+1 


P (cos0) 
p 


OO 


+ u2Jr 


2 £ 
P^2 


K P1(cos0) 

P P 

(p-2) (p+3) 


(D-19) 
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J I 


* r 


AIII  = AIV  at  r * R3 


ii&/A*iia,Avfc 
. — _ _ (rATT)  » — — (rA_)  at  r ■ R. 

P2  r dr  II  r dr  I 1 


1 i3/a  1 1 d , . x ^ 

7 57  (rAiii>  ■ s;  7 57  <rAn>  at  r - *2 


(E-2c) 


(E-2d) 


(E-2e) 


1 1 


r £ (rV  ar  7 Jr  (rAIII)  at  r = *3 


(E-2f ) 


After  appropriate  substitutions  of  equations  (E-la  to  E-ld)  into 
equations  (E-2a  to  E-2f)  the  following  boundary  value  equations  are  obtained. 


B 


A_ , R,  - A „R 


P+^2 


pl^l  P2  1 T u (p+1) 


(E-3a) 


Ap2R2P  + Bp2R2"('H'1)  " Ap3R2P  + “pS^'^1'  +(v2™2\  /p~2)  (p+»j 

(E-3b) 


S*P  * Rp3R3-<P+U  * ’ VR3-<Ptl) 


r2  [Ap2<R+1)Ri<p'1)  -<p)»p2Rr(p+2)]  - ^ APi(p+1^Ri 


(E-3c) 


(p-1) 


(E-3d) 
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I 


J 

J 


< 


II 


n 
! , 


“ i 


-(P+2) 


(E-3e) 


V2  [Ap3<P«>B2<P"1>  (p-2)(p^)  ] ' <E-3e) 

T2  [\2<',+1)R2<','1)  -<R>Bp2R2'<,>+2)] 

v2  [ -(R>VR3'(P+2’]-  T2  | (E-3f) 

+ ^3  P2JR3Kp  /(p-2)  (p+3)V 

These  algebraic  equations  provide  six  simultaneous  equations  with  six 
unknowns,  and  can  be  solved  for  the  coefficients  (A  t and  Bpl>  by  algebraic 
manipulation. 

Solving  algebraically  equations  (E-3a)  and  (E-3d)  results  in  the 
following  solution  for  A „ in  terms  of  B 

p2  p2 


(E-3f ) 


- UK 


where  [ X]  = -R^ 


(2p+l) 


l'*(-fe)S)/H) 


The  solution  for  A^  from  equation  (E-3a)  is: 


Apl  " Ap2  + Bp2Rl 


-(2p+l) 


(E-4) 


(E-5) 


Also,  solving  equation  (E-3c)  and  (E-3f)  algebraically  for  mathematical 
expressions  describing  B^  gives  the  expression  for  A^  in  terms  of  B^ 
and  known  quantities  after  equating  the  functions  for  B^  ■ B^  from  each 
equation. 
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»P3  -j-V  [(P-DR^ 


3y2JR2K 


2) (p+3) 


where 


+ [ X ] (p+DR^^K^1  -(p)R2"(p+2)K» 

|(-p)R2'(p+2)-  [z]  [x]  (p+1)R2p-1  + (p)R2~(p+2)[z] 

...  [■•(*%] 


y,JK  R. 
_ _£ 2 


-p+2 


(p-2) (2p+l) 


-(p+1) 


[Z]  1 


([X]  R2P  4 R2  <P+1) 


y2JR2K 


II 


tr-2)(&V+  V "2 

([X]  */  + R^1’) 


After  a numerical  value  for  is  calculated  on  the  computer 
problem  the  numerical  values  for  the  other  coefficients  can  be 
from  the  following  equations. 


Bp 2 ■ Bp ^ [z]  4-  Kp‘  , (see  equation  E-8) 


(E-9a) 


(E-9b) 


(E-9c) 


(E-9d) 


(E-9e) 


for  a specific 
obtained 

(E-lOa) 
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* Kp  (see  equation  E-6) 


Ap2  = [x]  8^2  (see  equation  E-4), 


Apl  » Ap2  + Bp2R^  ^ (see  equation  E-5) 


(p+3)K 


(E-lOb) 


(E- 10c) 


(E-lOd) 


B , = A „R„(2p+1)  + B „ + } p — (see  equation  E-7)  (E-lOe) 


pA  = Ap3K3  + p3  + (p-2) (p+3) 
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Y7.  J.  AKDAHAZY 
By  direction 


The  identity 


f 

•'r* 


ppi  (cos9)  Pp(cos9)  sin9d9 


2 (p+m) I . 

(2p+l)  (p-m) ! Vp 


(43) 

was  used  to  determine  K (equation  42).  Setting  J (6)  to  a constant  J 

Y V 

we  can  express  equation  (42)  as: 


K = 


p 2p (p+1) 


/ 

n 


(cosO)  sin9d6 


(44) 


Since  the  current  in 


in  this  problem  extends  from^-  ~ ajto  + aj, 
ssion  for  K (equation  44)  may  be  written  as 


see  Figure  4,  the  expression  for  Kp  (equation  44)  may  be  written 

7T 


K = 2P+1 

P 2p(p+l) 


I2 

“ - 


(cost)  sin9d9  + 


L t ~ Cl 


A p 


(cos9)sin9d0 

(45) 


By  noting  that 


p1(-cos?)  * (-l)P-1  P1(cos9) 


(46) 


it  follows  that  the  associated  Legendre  functions  Pp(cos9)  are  even  func- 
tions with  respect  to  cos9  when  p is  odd.  The  expression  for  K may  be 
simplified  to 


K = 4^4 


L 


P p(p+l)  J n „ ~p 


P^cosB)  sin9d6 


(47) 


2 - a 


when  p is  odd  (Kp  = 0 when  p is  even)  after  utilizing  the  symmetry  of  J(9) 
in  Figure  5.  By  changing  the  variable  in  equation  (47)  to  u ' = cos9,  the 


integral  for  K mav  be  written  as 
P 


K = - 2-P+1 

P p(p+l) 


- sin« 

/ Pp(u')dy' 

•'0 


(48) 


19 


J9  09  if  9™ 


Figure  1-B  - Infinitesimally  Thin  Current  Band 

Now  the  solution  for  ^ j (equation  B-5b)  can  he  reduced  to  that  for 
the  filamentary  coil  and  the  answer  compared  with  that  for  the  same  problem 
from  a standard  text  (Jackson,10  page  144,  equation  (5.46)). 

The  coefficient  JK^  for  the  expansion  of  the  current 


JA  9)  =E  JK  P1  (cos  0) 
V p P V 


(B-6) 


for  the  current  filament  is  obtained  from  the  equation  (see  equation  B-5d) : 


JK 


2p+l 

2p(p+l) 


/ 


cos(^)  7 

=■  6( cos0)P  (cos0)d(cos9) 
R,  P 

cos(O) 


(B-7) 


1-  I 2£±i \ pl 

R1  \ 2p(p+l)  / P 


(0) 
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where 


J(6)  * 6(eos9) 

R1 


The  special  value  of  the  associated  Legendre  function  is  calculated 
from  the  expression  (where  m = 1) 


P (0)  =j  (—0  . _ X' f (p+m,  even) 


(2p(^):  (ry5)  r 


(B-8) 


= 0 


> (p+m,  odd) 


■his  expression  may  be  derived  from  the  series  expansion  for  the 
Lvgendre  function  (see  Smythe)^  and  letting  cos0  = p1  approach  0, 

Substituting  the  expression  for  JK  (equation  B-7),  into  equation 
(B-5b)  results  in  the  expression  for  the  vector  potentials  in  the 
external  region  for  the  filamentary  coil. 


.11 


CO  P,IP  (cosB) 

V -I  £ 

p=l  (2p) (p+1) 


on 


(0) 


(B-9a) 


•>  4 7T 

Now  equation  (B-9a)  must  be  multiplied  by  — to  change  from  SI  units 

to  Gaussian  units,  and  p^  must  be  set  equal  tp  1 as  p^  is  the  free  space 

permeability  which  equals  1 in  the  Gaussian  v’stem. 

After  making  these  substitutions  we  have 


\>n 


00  I P1  (cos?) 

±1  £ -E 

* f.=1  2p(p+l) 


P (0) 
P 


(B-9b) 


N..w  taking  the  first  term  (p=l)  for  use  in  the  comparison  with  the  result 

• , , m . 

! n Jackson  we  havt 
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>11,  Lx 

(p-1) 


(cos9) 


(B-9c) 


Equation  (B-9c)  agrees  exactly  with  the  first  term  in  the  series  by 
Jackson, where  - a in  his  notation,  except  for  a minus  sign. 

Ulien  the  (-1)  phase  factor  in  the  definition  of  the  Associated  Legendre 
functions  P^(eos0)  is  taken  into  account,  the  two  expressions  are  the 
same.  The  general  expressions  for  the  Associated  Legendre  functions 
P™(cos0)  in  this  work  do  not  contain  the  (-l)m  phase  factor  used  by 

Jackson. 

The  | term  for  the  current  band  problem  also  reduces  exactly 

tp=l) 

to  Jackson's  term  for  n=0  (Jackson,  page  144,  equation  5.46)  for  the 
potential  inside  the  loop  when  the  (-1)  phase  factor  is  again  taken  into 
account . 
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